Decay of multiply charged vortices at nonzero temperatures 
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We study the instability of multiply charged vortices in the presence of thermal atoms and find 
various scenarios of splitting of such vortices. The onset of the decay of a vortex is always preceded 
by the increase of a number of thermal (uncondensed) atoms in the system and manifests itself by 
the sudden rise of the amplitude of the oscillations of the quadrupole moment. Our calculations 
show that the decay time gets shorter when the multiplicity of a vortex becomes higher. 



Experiments with atomic Bose-Einstein condensates 
have already shown their peculiar response to rotation, 
manifested through an induced irrotational flow and lead- 
ing to quantized vortices and superfluidity [TJ [2J |3J 0]. It 
was directly demonstrated by using the intcrfcrometric 
technique that the circulation is indeed quantized [5] and 
different aspects of dynamics of quantized vortices were 
investigated [H [HI [9] ■ Recently, also the persistent 
flow of condensed atoms in a toroidal trap was observed 

a- 

So far studies of quantized vortices were related to the 
single vortex state, the vortex lattices, and the multiply 
charged vortices. They included the analysis of the pro- 
cess of their nucleation as well as the stability conditions 
and the decay. For example, in Ref. [5] the lifetime of 
the single- vortex state in an axisymmetric trap is investi- 
gated at two different condensate parameters. Although 
in both cases the uncondensed part of the atomic cloud 
is almost undetectable its influence on the lifetime of the 
vortex is huge. Clearly, reducing the level of thermal 
atoms makes the lifetime of the vortex longer. Surpris- 
ingly long vortex lifetimes (up to several seconds) were 
observed for highly ordered triangular vortex lattices con- 
taining more than 100 vortices [3]. Multiply charged vor- 
tices have additional degree of freedom, their lifetime is 
limited by the instability leading to the decay into vor- 
tices having smaller charges. Pioneering work on doubly 
quantized vortices generated with the help of topological 
phase-imprinting technique showed that their lifetime is 
a monotonic function of the interaction strength [HUH]- 

Not much experimental work has been devoted to the 
rotational properties of Bose-Einstein condensates at fi- 
nite temperatures. In Ref. [7] the crystallization and 
the decay of vortex lattices is studied in the presence of 
thermal atoms. The decay of a lattice was observed in 
a nondestructive way just by monitoring the distortion 
of the rotating condensate and dramatic dependence on 
temperature was found. It turned out that the rotational 
frequency of the lattice decreases exponentially what has 
been attributed to the observation that the thermal cloud 
was also rotating. In contrast, when the thermal cloud 
is static the decay shows nonexponential behavior [2]. 
In Ref. [4 J , on the other hand, the influence of the ther- 
mal fraction on the stable circulation of condensed atoms 



was investigated. It was found that in a trap with ring 
geometry the persistent flow is possible even with the 
condensate fraction as small as 20%. 

Contrary to the case of zero temperature (for a review 
see |10j and references therein) there is only a small num- 
ber of theoretical papers considering the vortex dynamics 
at the presence of thermal atoms [11] . 

In this Letter, we focus on the evolution of multiply 
charged vortices generated in a way described in [1] and 
the role the thermal fraction plays in the splitting pro- 
cess. Therefore, we use an approach that allows a unified 
treatment of both condensed and thermal atoms. In this 
method, called the classical fields approximation [IS], the 
primary object is the complex field which represents all 
atoms. It is only the measurement process thereupon the 
atomic cloud is split in two components: the condensate 
and the thermal atoms. Such a decomposition is inher- 
ently built in the detection process since the detection 
always takes finite time and is performed under limited 
spatial resolution. In other words, the measurement is 
a kind of coarse graining procedure, hence eliminating 
some information from the classical field and allowing 
for the appearance of the thermal cloud. On a numerical 
level, it requires to calculate the time and/or space aver- 
age of a one-particle density matrix built of the classical 
field. Its dominantly populated eigenmode is the conden- 
sate wave function, the other modes represent thermal 
atoms, and its diagonal part integrated along the direc- 
tion of imaging beam is just what is monitored by the 
CCD camera. 

Our numerical procedure takes the following steps. 
First, we find the ground state of the Bose-Einstein con- 
densate in the toroidal trap made by combining the usual 
harmonic trap with the Gaussian laser beam (as in the 
experiment of Ref. [4]). The harmonic trap frequencies 
are taken as uj z = 2ir x 25 Hz and ui± — 2ir x 36 Hz. 
The blue-detuned laser beam serves as an optical plug 
repelling atoms from the trap center and the poten- 
tial from this beam is given by the Gaussian function 
Vo exp [— 2(x 2 + y 2 )/w ] (the laser beam propagates in z 
direction). Here, Vq is the maximum optical potential 
(Vo/h = 3600 Hz, where h is the Planck constant) and wo 
(u)q = 15 fj,m) is the waist of the beam. Next, we spec- 
ify the initial condition for the classical field by putting a 



particular amount of energy into the system and allowing 
the gas to thermalize. Since the classical field fulfills the 
time-dependent Gross-Pitaevskii equation [T^], we find 
its evolution and built the one-particle density matrix at 
each time. Now, we can do a coarse graining and calcu- 
late a fraction of condensed atoms, the condensate wave 
function, and the distribution of thermal (uncondensed) 
atoms. 

A recently developed experimental technique, based 
on a stimulated Raman process with Laguerre-Gaussian 
beams |H Q2] , allows to generate multiply charged vor- 
tices in atomic Bose-Einstein condensates. However, such 
vortices are unstable against the decay into a number of 
vortices with lower charges. It happens because: (i) the 
energy of a collection of vortices with lower charges is 
smaller than the energy of a single multiply quantized 
vortex, and (ii) the total circulation is preserved. Now, 
the question can be raised how the multiply charged vor- 
tex is split into a larger number of vortices. Since the 
energy of various configurations of vortices changes, one 
might expect different scenarios of a decay. 

To investigate the dynamics of the Bose gas after im- 
printing the phase on it we evolve the classical field ac- 
cording to the Gross-Pitaevskii equation. Before we go 
into details we want to stress that various scenarios of 
vortex decay are possible. The way the system follows 
depends on how fast the plug is switched off, i.e., on 
the energy of the gas when the plug is off. The shorter 
removal time ends in higher energy state, therefore al- 
lowing the access to higher energy configurations. For 
example, assuming the fivefold charged vortex was im- 
printed on the condensate, the considered configurations 
could be the one with four singly charged vortices placed 
on a circle around the vortex located in the center of the 
trap (called '4+1' configuration) or the other arrange- 
ment with five vortices settled on a circle and no vortices 
in the middle ('5+0' configuration). It can be checked nu- 
merically by using the imaginary time technique in a ro- 
tating frame of reference [Hj that the '4+1' structure has 
higher energy than '5+0' one when it rotates with large 
enough frequency (> 0.7 of radial trap frequency) and 
the opposite is true for slower rotation (even more, for 
slower rotation the obtained configurations are no longer 
at local energy minima). Therefore, depending on the 
time the plug is turned off the energy related arguments 
suggest possible scenarios of decaying process. Experi- 
mentalists could observe the '5+0' structure for slower 
removal of the plug and '4+1' configuration followed by 
the transition to '5+0' structure when the plug is taken 
off quickly. This scenario gets simplified if even as small 
as 8% fraction of uncondensed atoms is initially present 
in the system. In this case the system goes directly to 
'5+0' configuration. The first notice on the observation 
of different vortex configurations has been already given 

Our numerics confirms just discussed scenarios. Fig. 
[TJ\ shows the isodensity plots at times when the config- 
urations '4+1' (left frame, here the plug is off in 30 ms) 




FIG. 1: (color online). (A) Isodensity plots, corresponding 
to the density 1.2 x 10 13 cm -3 , for fast (30ms, left column) 
and slow (200 ms, right column) removal of the plug show- 
ing typical vortex arrays: '4+1' (left frame) and '5+0' (right 
frame) configurations. (B) The z-integrated total densities at 
times and conditions as in (A). (C),(D) The corresponding 
condensate and thermal densities, respectively. The size of 
each frame is 55 x 55 /im. 



and '5+0' (right frame, with the plug taken off in 200 ms) 
are present. The successive rows display the density as 
it is imaged by the CCD camera along the direction of 
axis of symmetry (B), the ^-integrated condensate den- 
sity obtained according to the classical fields approach 
(C), and the density of thermal cloud (D). Frames (C) 
and (D) clearly indicate that thermal atoms are located 
in vortices cores. An exception is the vortex placed at 
the trap center in '4+1' vortex array which is initially 
empty. However, immediately it is filled in with uncon- 
densed atoms it becomes unstable and moves away off 
the center forming the '5+0' configuration. 

Details of the dynamics of fivefold charged vortex are 
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given in Figs. [2] and [3J Fig. [2] shows the fraction of 
uncondensed atoms appearing in the system in the case 
of fast (30 ms) and slow (200 ms) change of the trapping 
potential. Fast removal of the inner plug means stronger 
disturbance of the gas, hence more effective production of 
thermal atoms. As it was already discussed in Ref. [15] . 
the level of the thermal noise strongly influences the life- 
time of the vortex. The larger number of thermal atoms 
the shorter lifetime of the vortex. The fivefold vortex 
decays approximately 1 s after the plug is taken out (ver- 
tical line (A) in Fig. [2]) and since the fast removal of the 
plug ends in a higher energy state, the higher energy vor- 
tex configurations are accessible, in our case it is '4+1' 
configuration. Afterwords, the singly charged vortices 
settled on a circle increase their separation from the vor- 
tex located at the center of the trap and simultaneously 
the frequency of the 4- vortex array rotation decreases. In 
other words, the 4-vortex array looses its energy and the 
lost energy goes to the thermal atoms. At some time (line 
(B) in Fig. [2| the frequency the 4-vortex array rotates 
approaches the critical frequency (~ 0.7 of radial trap fre- 
quency) and the system jumps to the '5+0' vortex config- 
uration increasing further the fraction of thermal atoms. 
Finally, both for fast and slow removal of the optical plug 
the system ends with '5+0' configuration rotating with 
approximately the same frequency The energy differ- 
ence between the initially imprinted the fivefold vortex 
and the final '5+0' vortex array is transformed to the 
thermal modes as can be verified by noticing that the 
thermal fraction in both cases remains approximately on 
the same level (see Fig. |2j. 
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FIG. 2: Condensate fraction as a function of time for fast 
(30 ms, black dots) and slow (200 ms, dark dots) removal of 
the plug. Vertical lines indicate: the times when the '4+1' 
(A) and '5+0' (B) configurations appear in the case of fast 
change of the trapping potential and the time when the '5+0' 
(C) structure is visible for slow removal of the plug. 

Our calculations show that the scenario of decay of 
a multiple vortex depends strongly on whether initially 
the Bose gas has a distinguished thermal fraction. In 
fact, the way the system behaves determines the upper 
limit for a number of uncondensed atoms. We found that 
already for the condensate fraction 92% the '4+1' con- 
figuration is not present irrespective of the duration the 
optical plug is taken off. Fig. [3] shows the evolution of 
the system for slow removal of the plug when the initial 



condensate fraction equals no = 0.92. In addition, we 
plot the time dependence of the quadrupole moment de- 
fined as J(x 2 — y 2 )\ip(r)\ 2 d 3 r, where ^(r) is the classical 
field. It turns out that the behavior of the quadrupole 
moment exhibits clearly the time the multiple vortex be- 
gins to split . The decay of the vortex is accompanied by 
the sudden increase of the amplitude of the quadrupole 
moment and allows unambiguously to determine the life- 
time of the vortex for any vortex charge (see Fig. EJ). 
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FIG. 3: (color online). Decay of a fivefold charged vortex. 
The upper panel shows the density at various times (from left 
to right): while the plug is on (A), after the plug is off (which 
takes 200 ms) but before the vortex is split into the singly 
charged vortices (B), and when the vortex is clearly split into 
five vortices (C). In the lower frame the quadrupole moment 
is plotted as a function of time showing clearly the onset (0.3 s 
after the plug is taken off) of the decay process. 




FIG. 4: Quadrupole moment as a function of time for dif- 
ferently charged vortices (from doubly charged to sevenfold 
charged as indicated in the left upper corner of each frame). 
Vertical lines show the beginning and the end (delayed by 
200 ms) of removal of the plug. Note that, in fact, the dra- 
matic change in the behavior of the quadrupole moment de- 
termines the onset of the decay of a vortex. 

Filially, in Fig. 5 we plot the decay time of the multi- 
ple charged vortex as a function of its charge, assuming 
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FIG. 5: Decay time as a function of the charge of a vortex 
for a given initial fraction of condensed atoms (no = 0.92) for 
slow (200 ms) and fast (30 ms, Inset) removal of the plug. 
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of fast change (30 ms). First observation is that faster re- 
moval of the plug leads to shorter decay times. This fea- 
ture can be understood based on Fig. [2] Faster removal 
produces more thermal atoms and hence the destabiliza- 
tion of the vortex happens in a shorter time (see also Ref . 
[16]). The same arguments help to understand the de- 
pendence of the lifetime on the charge of the vortex. As 
it is displayed in Fig. [6] the production of uncondensed 
atoms is enhanced when the topological charge of the 
vortex is bigger. 

In summary, we have studied the decay process of mul- 
tiply charged vortices. We show that various scenarios of 
vortex splitting are possible depending on the level of un- 
condensed atoms appearing in the system as a result of a 
change of the trapping potential. Initially, an amount of 
uncondensed atoms is determined by how fast the plug 
supporting the toroidal trap is removed. Faster removal 
produces more thermal atoms and leads to quicker decay, 
however, afterwords possibly displaying different vortex 
configurations. Finally, both fast and slow removal re- 
sults in the same vortex configuration. At the time the 
vortex decays the large oscillations of the quadrupolc mo- 
ment appear. The decay time depends also on the mul- 
tiplicity of the vortex and for even a tiny uncondensed 
fraction gets shorter for higher multiplicity. 



FIG. 6: Condensate fraction as a function of time for fast 
(30 ms) removal of the optical plug. Squares, triangles, and 
diamonds correspond to the charge of the imprinted vortex 
equal to 5, 6, and 7, respectively. 

it is the time that elapsed since the plug beam was off. 
Here, the initial condensate fraction equals no = 0.92. 
The main frame corresponds to slow change of the trap- 
ping potential (200 ms), whereas the inset shows the case 
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